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Abstract 

We summarize the present state of research on the darkon fluid as a model for the dark 
sector of the Universe. Nonrelativistic massless particles are introduced as a realization of the 
Galilei group in an enlarged phase space. The additional degrees of freedom allow for a non- 
standard, minimal coupling to gravity respecting Einstein's equivalence principle. Extended to 
a self-gravitating fluid the Poisson equation for the gravitational potential contains a dynam- 
ically generated effective gravitational mass density of either sign. The equations of motion 
(EOMs) contain no free parameters and are invariant w.r.t. Milne gauge transformations. 
Fixing the gauge eliminates the unphysical degrees of freedom. The resulting Lagrangian pos- 
sesses no free particle limit. The particles it describes, darkons, exist only as fluid particles of 
a self-gravitating fluid. This darkon fluid realizes the zero-mass Galilean algebra extended by 
dilations with dynamical exponent z — | . We reduce the EOMs to Friedmann-like equations, 
derive conserved quantities and a unique Hamiltonian dynamics by implementing dilation 
symmetry. By the Casimir of the Poisson-bracket (PB)-algebra we foliate the phase space and 
construct a Lagrangian in reduced phase space. We solve the Friedmann-like equations with 
the transition redshift and the value of the Casimir as integration constants. We obtain a 
deceleration phase for the early Universe and an acceleration phase for the late Universe in 
agreement with observations. Steady state equations in the spherically symmetric case may 
model a galactic halo. Numerical solutions of a nonlinear differential equation for the grav- 
itational potential lead to predictions for the dark matter (DM) part of the rotation curves 
(RCs) of galaxies in qualitative agreement with observational data. We also present a general 
covariant generalization of the model. 

'Accepted for publication by 'Entropy' in the special issue 'Modified gravity: From black holes entropy to current 
cosmology' 
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1 Introduction 



Understanding the dark sector of the Universe is one of the greatest challenges of today's 
theoretical physics. The main question is whether the underlying astrophysical observations 
(missing gravitation in galactic systems resp. the accelerated expansion of the Universe) 
can be understood by means of known physical concepts, e.g. General Relativity (GR), or 
whether some kind of "new" physics is needed. Indeed, cosmic acceleration could be an 
apparent effect due to the averaging over large scale inhomogeneities in the Universe, but a 
reliable quantification of such an effect is not yet available (see the recent review [I]). On 
the other hand there exists an overwhelming evidence for the existence of gravitational effects 
on all cosmological scales (termed "dark matter" (DM)) which cannot be explained by the 
gravitation of standard matter in the framework of GR (see the review [2]). 

The current picture of cosmological structure formation assumes for DM some pressure less 
dust containing massive, perhaps weakly interacting cold particles (CDM). The present stan- 
dard cosmological model (ACDM-model) uses besides CDM a positive cosmological constant 
A as the cause for the accelerated expansion of the Universe (see recent reviews on dark energy 
(DE), e.g. [3J, [J and [S]). But, "favoured by a number of observations" [3j, the ACDM-model 
suffers at least from the following insufficiencies: 

• Interpreted as the energy density of the vacuum, the experimental value of A turns out 
to be a factor of about 10 54 too small (see the recent estimate in [B]). 

• None of the proposed CDM-constituents has been observed (cp. [?])• 

• Rotation curves (RCs) of DM-dominated galaxies behave in their inner part in sharp 
contrast to the CDM-based simulations (known as the core-cusp problem [8]). 

In a very recent paper P. Kroupa [9] lists a large set of extragalactic observations which 
falsify the ACDM-model. 

In this paper we will neither comment on the huge number of phenomenological models 
for a dynamical DE (see [3]- [5]), which all rest on at least one unknown function {e.g. a scalar 
field potential), nor on those modified gravity theories which explain DM-effects only [e.g. 
modified Newtonian dynamics (MOND) in its simplest, Milgromian form [IQ]). But we know 
of at least two kinds of modified gravity models which explain both, DE and DM i.e. the 
whole dark sector of the Universe: 

• For MOND-models see section 9 of the very recent review [TT]. But all these models are 
very phenomenological in that they are based on one unknown "MOND function" . 

• Modified theory of gravity (MOG) by Moffat |12| in which, besides the Einstein-Hilbert 
and matter actions, a massive vector field is introduced whose mass, coupling constant to 
matter and gravitational constant are promoted to scalar fields. The self-interaction po- 
tentials of these four fields are not fixed a priori and play the role of DE [13] . This model 
can explain galactic phenomena as e.g. RCs Q3] as well as the accelerated expansion of 
the Universe and other cosmological observations 13 . 

In summary, all present attempts to explain quantitatively the new cosmological obser- 
vations [e.g. late time cosmic acceleration) as well as the galactic RCs and other galactic 
phenomena contain either some new parameters or even free functions. If we assume that 
the cosmic acceleration is a real effect and not an apparent one (see above) then, obviously, 
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we need some new (i.e. unconventional) physics which, however, should be based on known 
fundamental physical principles (e.g. symmetry). 

In this paper we review some ongoing research on a model which is a first building block 
for a new theory describing the dark sector of the Universe. This model, introduced in [15] and 
further developed in [T§] and Q~7], contains no new parameters in its Lagrangian and is based 
on well- known physical principles (Galilei symmetry and Einstein's equivalence principle). 
The only free parameters appearing in our model are some integration constants (see sections 
6 and 7). The basic idea it uses involves nonrelativistic massless particles. The use of such 
particles may appear very strange at a first sight as it seems to contradict Special Relativity 
according to which massless particles must move with the velocity of light. However, as 
we have shown before, we can introduce nonrelativistic massless particles if we enlarge the 
dimension of phase space; i.e. such particles are not described by position and momentum 
only. The additional degrees of freedom allow for a nonstandard, minimal coupling of these 
"exotic" particles [15] to gravity in accordance with Einstein's equivalence principle. This 
coupling, for a self-gravitating fluid, leads to a dynamically generated effective gravitational 
mass density of either sign which is the source of the corresponding gravitational field. This 
fact leads to the possibility of using such a model as a first building block of a new theory 
describing the dark sector. In our approach we formulate the equivalence principle (local 
equivalence between gravitation and acceleration) as usual as an invariance of the dynamics 
w.r.t. arbitrary time-dependent translations (Milne gauge transformations). By fixing the 
gauge we can reduce the dimension of phase space and remove the unphysical variables. The 
Lagrangian of the resultant model does not possess a free particle limit and hence the particles 
it describes, called darkons, exist only as fluid particles of a self-gravitating fluid. This darkon 
fluid presents a dynamical realization of the zero-mass Galilean algebra extended by anisotropic 
dilations with the dynamical exponent z = |. This dilation symmetry emerges in a natural 
way from the minimal gravitational coupling. The cosmological EOMs are derived, as usual, 
from the darkon fluid EOMs by restricting the form of their solutions through the cosmological 
principle. However, the corresponding Hamiltonian cannot be derived from the darkon fluid 
Hamiltonian. So we have to derive a Hamiltonian and the corresponding Poisson-brackets 
(PBs) from the EOMs, a procedure which, as is well known, is not unique. So we construct, 
using the constants of motion, a Hamiltonian which has the correct scaling dimension and 
obtain the unique PBs by implementing the dilation symmetry. Moreover, our procedure 
allows us to construct a Lagrangian having foliated the phase space by means of the Casimir 
of the PB-aLgebra. We solve the cosmological EOMs and obtain, by fixing two integration 
constants by their values determined from experimental data, a prediction for the Hubble 
parameter H(z) (z is the redshift) in general agreement with the existing data. 

Furthermore, we demonstrate that our model not only gives a proper description of the 
late time cosmic acceleration but it also predicts galactic halos and, qualitatively correct, the 
corresponding flat rotation curves (RCs). To do this we derive for the spherically symmetric 
steady state case a nonlinear ordinary differential equation for the gravitational potential 
which is then solved numerically. 

Almost all cosmological models are based on General Relativity (GR) or on some of its 
modifications which, however, preserve general covariance. The appearance of anisotropic 
dilations in our model seems, at first sight, to forbid a generalization of our model to a general 
covariant theory. But the dilation symmetry with z = | is an emergent symmetry and not an 
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input of our model. So we are able to propose a generalization of our model showing general 
covariance without imposing any form of scale symmetry. 

One of the aims of the present paper is to show in detail the main ideas of our model. A 
comparison with observational data (Hubble parameter resp. galactic rotation curves) serves 
only as an illustration that it might be meaningful to develop the present model to a proper 
theory of gravitation for the dark sector of the Universe. For this reason we do not produce 
any least square fits for the undetermined integration constants appearing in the cosmological 
context . 

The paper is organized as follows: In section 2 we introduce Galilean massless particles and 
their coupling to gravity. In section 3 we describe the generalization of this particle picture 
to a self-gravitating fluid, fix the gauge and so obtain the darkon fluid. Then, instead of 
the Poisson equation, we get a transport equation for the gravitational field. In section 4 
we describe the Hamiltonian dynamics of the darkon fluid to have a basis for quantization. 
Space-time symmetries of our model are discussed in section 5. The darkon fluid cosmology 
is presented in section 6. Starting with the cosmological principle we derive Friedmann-likc 
equations, discuss the corresponding Hamiltonian dynamics and compare our predictions with 
the existing data for the Hubble diagram. In section 7 we consider steady state equations in the 
spherically symmetric case and derive a nonlinear, ordinary differential equation (ODE) for the 
gravitational potential. Numerical solutions of this ODE lead to predictions for the behaviour 
of the DM-part of galactic RCs. In section 8 we outline how to determine a correction to 
Newton's gravitational law due to the cosmic expansion. Finally, in section 9 we present a 
general covariant generalisation of our model. We close with some conclusions and give an 
outlook for future work. 

2 Galilean massless particles and their nonstandard 
coupling to gravity. 

As stated earlier the introduction of nonrelativistic massless particles may appear very strange 
at first sight. It seems to contradict Special Relativity according to which massless particles 
must move with the velocity of light. But in Special Relativity spinless particles are described 
by the vectors of position and momentum only. In our recent papers |15j-|17] we have shown 
that we can introduce nonrelativistic massless particles if we enlarge the dimension of the phase 
space from six to twelve. In subsection 2.1 we will show that this enlargement of phase space 
is also necessary: The 12-dimensional phase space is the minimal one allowing the dynamical 
realization of nonrelativistic massless particles. 

How to couple gravity if there is no mass? Our nonrelativistic framework forbids the 
coupling to energy, as then this would be a relativistic effect. So we take Einstein's equivalence 
principle as a starting point in subsection 2.2. Then, within the Lagrangian framework, the 
gravitational field will be coupled to one of the new degrees of freedom. 

2.1 Galilean massless particles 

Our first task involves finding a minimal dynamical realization of the unextended (m = 0) 
Galilean algebra Go in three space dimensions. 



5 



Go is given by the Lie-brackets 



[A iy Lj] = UjhAh, At G (Pi,Ki,Li), 

[H,Ki]=Pi, [H, Pi] = [H, Li] = 0, (1) 

[Pi, Kj] = [Pt, Pj] = 0, [Ki, Kj] = 0, 

where H,Pi,Ki and Li are the generators of time-resp. space translations, boosts and rota- 
tions. 

A minimal dynamical realization of |T]) involves a construction of the elements of Go in 
terms of a minimal number of phase space variables and without the introduction of any free 
functions or parameters. Then the Lie-brackets in |T]) have to be understood as Poisson- 
brackets (PBs) w.r.t. the PB-algebra of phase space variables. 

Let x(t) be a particle trajectory and pthe corresponding conjugate momentum satisfying 
the fundamental PB-relation 

[Xi, pj] = 8ij (2) 
The translation generator P acts on x resp. p as 

[xi, Pj] = Sij resp [pi, Pj] = 0. 

Therefore we get the identification 

Pi = Pi. (3) 

To show that our minimally enlarged phase space contains, besides x and p, two additional 
vectors corresponding to the reduced boost q and the velocity y, we proceed in two steps: 

• 1. Define a reduced boost vector q 

qi = Ki - Pit. (4) 

The boost K acts on x as 

[ Xi , Kj] = Sijt (5) 

leading by Q and Q to 

[an, qj] = 0. (6) 

Taking the PB of Q with Pj it is easily seen by the expressions in |T]) that the PB of Pi 
and qj vanishes 

\pi, Qj] = 0. (7) 
Despite of (J©J) and |7| q turns out to be nontrivial: Taking the PB of Q with H gives 

<ji = [qi, H] = - Pi . (8) 

Thus we conclude from |6|-([8]l that the qi are independent dynamical variables which, 
according to Q and Q, are commutative 

[qi, ©] = 0. (9) 



(i 



2. Next we introduce the velocity vector 

= [xi, H], (10) 

which is translationally invariant (take the PB of ( |10[ ) with pj and use the Jacobi identity) 
to get 

[Vi, Pi] = 0. (11) 
From flOjl, 0, j§ and jHJ we get the action of the boost generator on 

[W, «i] = (12) 
Inserting now (|4| into ( 12 1 leads by (111 to 

[Vi, «j] = (13) 

We conclude from ( 13 \ and by |[6]|, Q and Q that ?7 cannot be a function of the other 
variables x, p or <f.) 

To close the PB-algebra in phase space we choose, in accordance with the principle of a 
minimal dimension, 

[ Xi , Vj ] = 0. (14) 



Note that ( 14 1 respects all Jacobi identities. 

It remains to construct out of the phase space variables x, p, q, and y the Hamiltonian H 
and the generator of rotations L. 

Our foregoing results lead to the EOMs 



x% = Vi, Pi = and g 4 = -p t 
but iji has not not yet been fixed. 



(15) 



The EOMs ([15} together with the canonical PBs 

[xi,pj] = Sij, [yi, q-j] = 5ij (16) 
imply that the Hamiltonian is given by 

H = piyi + f(q) 

The principle of a minimal dynamical realization of Go forces us to put / = 0. Thus we see 
that the Hamiltonian H is given by 

H = VlVl , (17) 



which leads, due to (16 1, to the EOM for yi, namely 

m = o. 



(18) 



The conserved generator of rotations Li is then given by 

L t = dihi (x k pi + ykqi)- (19) 
The two parts of Li act separately on the (x,p), resp. (y, q) parts of the phase space. 



The Hamiltonian dynamics given by H (17 'I and the nonvanishing PBs (161 is equivalent 
to the dynamics derived from the Lagrangian 

io = Pi±i + qiiji — piyi. (20) 

Finally we conclude 
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• A minimal dynamical realization of the unextended Galilei algebra in three dimensional 
configuration space requires a 12-dimensional phase space spanned by (x, p, q, y) 

• The relation between energy H and the momenta (velocities) is nonstandard, 

• Momenta pi and velocities yi are dynamically independent, 

• The qi have no counterpart in standard particle mechanics, 

• Galilean massless particles move with arbitrary finite velocity. 
Dilations 

Go may be enlarged by dilations with arbitrary dynamical exponent z. The corresponding 
generator has to satisfy the following PB-relations (cp [18] for z = ^r) 

[D,H] = -H, [D, Pi] — ~^Pi> (21) 

[D,Ki] = (l-^jKi, [D, Li] = 0. 
It can easily be checked that D is given by 

D =tH - ^x iPi + fl - i j y iqi . (22) 

Below we will show that the value of z becomes fixed by the coupling of our massless 
particles to gravity. 

2.2 Coupling to gravity 

To couple our massless particles to the gravitational field strength gi{x,t) we have to be 
consistent with Einstein's equivalence principle: Locally, the gravitational field strength is 
equivalent to an accelerating frame. The only known nonrelativistic EOM for the particle 
trajectory satisfying this form of the equivalence principle is given by the Newton law: 

Xi(t) = gi (S(t),t) (23) 



because (231 is invariant w.r.t. arbitrary time-dependent translations (Milne gauge transfor- 
mations) [19] 

Xi ->■ Xi = Xi + a,i(t) (24) 

provided that gi transforms to 

g'i{x,t) = g l {x,t) + ai(t). (25) 



The EOM ( 23 \ is realized if we add to Lq an interaction part (minimal coupling) 

L tnt = ~qigi(x,t) (26) 
Then the EOM q\ — gets replaced by 

q\ = qkdigk- (27) 
The total Lagrangian L = Lo+Li nt now becomes invariant w.r.t. the gauge transformations 
Xi = Xi + cn(t), ,Ui = yi + a,i(t), Pi = pi, qi = qi, (28) 
9i(x ,t) = gi{x,t) + oi(t). 
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This gauge invariance, resp. Einstein's equivalence principle still holds if we add in ( 26 1 
further terms containing spatial derivatives of gt i.e. by replacing 

gi — > g% — KiAgi + higher order terms, (29) 

where A is the laplacian and K\ is a constant. 

3 A self-gravitating (darkon) fluid 
3.1 Lagrangian formulation 

In the last subsection (2.2) we have introduced a model for the motion of one Galilean massless 
particle in an external gravitational field g(x,t). 

Now we want to generalize this model in two respects: 

• 1. Instead of one particle we consider a continuum (fluid) labeled by comoving coordi- 
nates £ £ R 3 (Lagrangian formulation of fluid dynamics, cp. [2U] ) . 

• 2. We treat the gravitational field g as a dynamical variable in order to get a self- 
gravitating fluid. 



To realize the first point we substitute in the EOMs ( 23 1 and ( 27 1 



Xi{t) -> Xi{£,t) and qi (t) ->■ qi{i,t). (30) 

To realize the second point we add to our previous Lagrangian Lo + Li nt a kinetic part 
L field for the gravitational field. We choose as usual 



1_ 

Afield - ^ 



d 3 x( 9i (x,t)) 2 , (31) 



where G is Newton's gravitational constant. 

Then our Lagrangian, formulated in terms of (xi,Xi), {qi,qi) and g, becomes 



L = 



-n J d 3 £ t)±i(£*) + qi(Ct) gi (x(lt),t)) + Lf ie id, (32) 



where no is the (constant) particle density in £-space. 

Varying q; , resp., Xi leads to the continuum generalization of the previous EOMs (231 and 
27b 



Xi{i,t) = gi (x(Z,t),t) (33) 

and 

ft (6 1) = q k & t ) A g k {x(C t ) , t ) . (34) 

Now we have two possibilities to go further: 

In the first case we take gi = —di(j> and vary <j> getting 

= 4TTGd k (nq k ), (35) 

with 

n(x,t) = no J _1 ||*_|*(g t ) (particle density) (36) 



and 



J = det I dX l^ f) I (.iMc.bian). C!7) 

Ot,k 
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where £(x,t) denotes the function inverse to af(£,t). 



The right hand side of the Poisson equation (35 1 describes a dynamically generated effective 
gravitational mass density which may be of either sign with 
+ sign leading to an attractive gravitation, 
- sign leading to a repulsive gravitation. 

This promotes the self-gravitating fluid to being a possible candidate for the dark sector 
of the Universe. 



Note that our EOMs are invariant w.r.t the gauge transformations (see (241, (25l) 

Xi -> Xi + a,i(t), (38) 

when 4> transforms as 

4>(x,t) -» 4>'{x' ,t) = <j){x,t) - &i(t)xi, (39) 

and qi and n remain invariant. 

In the second case we vary <?j instead of <j> and this leads to a linear relation between qi 
and 

g l (x,t) = -4tt Gn(x,t) qi(£{$,t),t). (40) 
Note that the expression for gi can be rewritten as 



(41) 



The expression (41 1 allows to eliminate, from the EOMs, the unphysical phase space vari- 
ables qi in favour of the gravitational field strength gi. 

This removes the unphysical degrees of freedom. However, we have to pay a prize for this 
as now 

• the gravitational field strength is not automatically given by the gradient of a potential, 

• the expression (41) destroys the gauge symmetry i.e. it fixes the gauge. 
Remark The gauge symmetry can be restored by considering a curl-free velocity field 

u — Vu 

and treating u as the variable to be varied |16j . 

Inserting qi from (41 1 into the previous EOMs we get two coupled EOMs for x(£,t) and 

g{x{lt),t) 



dt 2 



Xi{£,t) = gi{x{£,t),t) 



(initial condition x(£, 0) = £) and 



dt 2 



(I) - s** 



Furthermore by inserting (411 into the Lagrangian (321 we obtain 

1 



L = 



no f 3t d / ffi (£(£,*),*) 



4ttG 



d^ 



dt 



Xi(£,t) + 



8ttG 



d xgi (x,t), 



(42) 



(43) 



(44) 



where we have used the identity 



d 3 £Jg 2 i{x{£,t),t) = / d 3 xgi(x,t) 
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It can be easily checked that the Euler-Lagrange EOMs following from (44 1 are identical 



with (421 and (431 



The EOMs do not contain the gravitational constant G. The Lagrangian contains it only 



as a common factor. Moreover, our modified Lagrangian (44 1 does not split into two parts 



involving a free and an interacting term. Clearly this strange property is a consequence of 



the phase space reduction due to (41 1. We know only one other physical system possessing a 



similar property: the interacting Chaplygin gas if one eliminates the particle density from the 
Lagrangian (see [21] . Sect. 2.1, item (i)). 

3.2 Relabeling symmetry and the transport equation for the 
gravitational field 

Any fluid dynamics in the Lagrangian formulation has to be invariant under infinitesimal 
relabeling transformations (volume preserving diffeomorphisms) (see |20| ) 



£ + a(|), with V • a = 0. 



For the corresponding Noether charge Q we obtain from the Lagrangian ( 44 \ 



fffc dx k 
n d^i 



(45) 



(46) 



where we have alerady performed one partial integration. But from the EOMs (421, (43 1 it 



follows that the integrand in 1 46 1 is already conserved i. e. that 



dx k 5* 
d£i n 



dx k /pfe 
d£i dt V n 



(47) 



is a conserved quantity and so is a function of £ only. By solving (471 for (— ) we obtain 



the once-integrated form of the EOM ( 43 1 

(?)- 



d_ 

dt 



Qk dx k 
n dxi 



0*(O 



dxi '" 



(48) 



which is a transport equation for the gravitational field (cf. [22] for the standard self- 
gravitating dust fluid). 



3.3 Eulerian formulation 

In the Lagrangian formulation the fluid particles are labeled by comoving coordinates £. So 
a Lagrangian observer moves with the particle. On the other hand an Eulerian observer is 
located at a fixed position x (see [22]). Therefore we obtain the fluid fields in the Eulerian 
formulation from the ones given in the Lagrangian formulation by the following change of 
independent variables 

(£,t) -> (x,t) with x = x(i,t). (49) 
Correspondingly, we have to transform the time derivatives 

^- — > Dt — tt + u k (x, t) -J^— (convective derivative) (50) 
dt at ox k 

with 

u k (x,t) = — x k (i,t)\^ 3)ty (51) 

The EOMs then become: 
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Euler eq. 

D t Ui = gi. (52) 
Transport eq. for the gravitational field: 

D t (^) = 9 -3 % u k - 9i, (53) 

where the auxiliary variable 
obeys the EOM 

D t 6i + 6 k di u k = (54) 
and the particle density n obeys the continuity eq.: 

d t n + d k (nu k ) = 0. (55) 

All these equations can be derived from the Lagrangian 

L = J ' d 3x (gi(D t Ui - ^gi) + 6(d t n + d k (nu k )) - naDt/3^ . (56) 

Here 9, a and /3 give the Clebsch-parameterization of the vector field 9i(x,t) 

9 % = d l 9 + adiP. (57) 

They obey the Euler-Lagrange EOMs 

D t = D t a = D t /3 = 0. (58) 



Operating with D t on (57 1 leads then to the EOM (541 for di. All the remaining EOMs are 



Euler-Lagrange EOMs derived from ( |56[ |. 

3.4 Coupling with baryonic matter 

We describe baryonic matter (B) by a standard self-gravitating fluid coupled to the darkon 
fluid by the jointly created gravitational field (j>- Then, within the Lagrangian formulation, we 



have to add to (32 1 the following baryonic Lagrangian 



Lb = m J d 3 e ( (±f( ^ ))2 - ^ B (£t)) ! . (.I!)) 

where m is the (constant) baryonic mass density in £ -space. 
Varying xf leads to the Newton equation 

xf(£,t) = -d i <t>(3S B (£,t),t), (60) 

whereas varying <f> in the total Lagrangian 

L — L D + L B + Lfi e id, (61) 



where Lr> is the darkon fluid Lagrangian (first term in ( 34 1 ) , leads to the Poisson equation 



V<j> = 4ttG (d k (nq k ) + p B ) (62) 

whose r.h.s. consists of the sum of the darkon fluid source and the usual _B-source. The B-mass 
density pB is defined by 

p B (x,t) = m J _1 |f =fet) , (63) 



12 



where J is the Jacobian defined by (37 1. We note that ps satisfies the local conservation law 
(continuity equation) 



d t pB + d k (u k p B ) = 



(64) 



where u a is defined analogously to 11511. 

Adding baryonic matter leads to the following changes in the foregoing subsections: The 
gravitational force appearing in Lb is the total one g$ = —di<f> and so is it in the case of 



the EOMs (33 1 and (34 1. The solution of the Poisson equation (621 is a superposition of the 
gravitational potentials created by the B and the D-sources 



with 



and 



Ac 



A-kGpe 



dk <t>D = A-KGnq k . 



(65) 



(66) 



(67) 



So (|67|) is now the substitute for (41 1 and, therefore, the conserved quantity (47 1 becomes 

>eling transfor: 

dx k d ( d k <t>D 



(note that Lb is separately invariant w.r.t. relabeling transformations ( 45 1 ) 

. dx k d k <t>D 



+ 



9& dt 



(68) 



Then, in the Eulerian formulation, we obtain finally instead of ( 53 1 the transport equation 

' di4>D\ d k (f>D 



d t u k + Oi 



(69) 



where (9; obeys the EOM ((54 



4 Hamiltonian dynamics 



The derivation of the Hamiltonian dynamics, le the formulation of classical dynamics in terms 
of a Hamiltonian and a PB-algebra, is an indispensable prerequisite for any quantization pro- 



cedure. We start with the Lagrangian ( 56 1 from which we obtain immediately the Hamiltonian 
H = J d 3 x ^gi(u k d k Ui - igi) - nO % u^ , (70) 

where we have used the identity [21] 

9d k {nu k ) - nu k ad k p = -n9 k u k + d k (u k n0) (71) 



and then neglected the boundary contribution arising from the last term in (71 1. For the 
Hamiltonian EOM 

A(x,t) = [A(x,t), H] (72) 



we choose as independent field variables A £ (n, 0, net, ft, V4, gi). To read off from (701 and 



( 72 1 the PBs for the A-fields we next rewrite their EOM as follows 

n = -d k (u k n), 6 = -u k d k 6, 
(na)' = -d k (u k na), $ = -u k d k (3, 



Ui — —Uk d k Ui + gi and (73) 
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where 



9i = gkdiUk - d k (u k gi) - n6i 
9, = d % 9 + adi/3. 



The EOM ( 73 I have the Hamiltonian form ( 72 1 if the fluid fields A obey the following 



nonvanishing PBs 



and 



[9{x,t), n(x ,£)] = AtyGS(x - x ) 
\P(x,t), (na)(x ,t)] = iwG6(x-x') 



(74) 
(75) 



[gi(x,t), uj(x ,t)} = 4n G 6 i:i 8(x - x ). (76) 

For these PBs all Jacobi identities are trivially fulfilled. It is interesting to note that the 
PB-algebra ( 74|76 1 decomposes into three irreducible parts. Each component has a canonical 
structure which allows quantization by applying Dirac's rule (details will be given in a separate 
paper) 



[A, B]i 



ih 



\A, B] 



(77) 



where A(B) denotes the quantum operator corresponding to the classical function A(B) . The 
PB-algebra ( |74|76 \ possesses a Casimir (being a quantity whose PBs with all elements of the 
PB-algebra vanish) defined by 

d6i(x, t) 



C 



d X£j, 



dx k 



■6i(x,t). 



To prove that C is a Casimir we proceed in two steps: 



We choose a particular relabeling transformation (451 with Oi(£) = — ei 



pom 

1 8£ k ■ 



corresponding Noether charge ( 46 1 is then given by 



(78) 



The 



(79) 



But Eulerian variables do not respond to a relabeling of the Lagrange variable £ [21) . So 
Q is a Casimir of the PB-algebra ( 74)76 1. 

• It remains to express Q in terms of Eulerian variables. Do do this we get from the 
definition of 9i(x,t) 



HO = t),t) 



dxk 



However, using ( 80 1 and the identity 

dXn dXm 



we obtain 



_ r % 

o^i d£i dxj 



Cjfcl^T — J tmjn ^ t>m \ X > t )\g=3(t,t) 



(80) 



(81) 



d£ k 



dxj 



Hence we finally have 



Q = C. 



(82) 



Of course the Casimir property of C can be shown also directly by bracketing ( 78 1 with 
all fluid fields and using the PB-algebra |74|76| . 



We recall that C vanishes if 9i(x,t) is given by the Clebsch-parametrization ( |57[ ) with 
non-singular functions 9, a, (3 vanishing at infinity [21j . 
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We note that the vortex helicity 

H e = J d 3 xu(x,t) ■ u(x,t), (83) 

where 

uj = V X u 

which is a Casimir in standard fluid dynamics, is not a Casimir in the case of our model. H e 
is a conserved quantity if g is the gradient of a potential, but it has a non- trivial PB with g, 
namely: 

\gi(x,t), H e ] = torGu)i(x,t). (84) 
Up to now we have used, within the Hamiltonian framework, those fluid fields which lead 
to the simplest possible form ( 74|76 I for the PB-algebra. But in many applications it is more 



convenient to introduce instead of the auxiliary field 6i the momentum field Pi 

47vG P i(x,t) = D t i^) = ^cU fe - 9 i; (85) 

where we have used the definition of pi in the Lagrangian formulation 

1 d ( gi s 



2(C,i) 



and then changed to the Eulerian formulation. In terms of the momentum density npi we 
obtain for our PB-algebra, besides the PB-relation (1761), 



[(npi)(x,t), n(x ,t)] = —n{x,t)diS{x — x), (86) 

[(npi)(x,t), u k {x ,t)] = diu k 5(x - x ), (87) 

[(npi)(x,t), g k (x ,t)] = -gk(x,t)di8(x — x), (88) 

[(npi)(x,t), {npj){x ,t)] = -(npj)(x,t)di8(x - x ) + (npi)(x ,t) d'j 8(x - x ), (89) 
which can be obtained from the transformation formulae 

n(x,t) = / d'^S(x- x(€,t)) (90) 



(nA)(x,t) = I d 3 £A(g,t)5(x-x(£,t)), Ae( Pl , Ul , gi ) (91) 



and the canonical PBs in the Lagrange formulation (continuum generalization of (16l) 

[xi(l t), Pi(t,t)] = M(f~ t) (92) 

and 

Mi,t), qj (t,t)] = M(f-f") (93) 

if, in addition, we use the relation ( |4"T| for gj . 

5 Space-time symmetries 

In this section we will show 

• The darkon fluid possesses the unextended (zero mass) Galilei algebra Go as a symmetry 
algebra. The corresponding energy-momentum tensor is not symmetric and contains a 
negative pressure term. 

• Go may be enlarged by dilations with dynamical exponent z = |. This particular value 



of z is determined by the minimal coupling to gravity as given by eq. ( 26 1 
Expansions are not a symmetry of the darkon fluid. 
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5.1 Space-time translations and the energy-momentum tensor 

The linear momentum (generator of space translations) is given by 

Pi = J d z x{n Pl )(x,t), (94) 

where the momentum density npi obeys the local conservation law 

1 2 

d t {npi) = -d k (u k npi) + —^,dig k = -d k T ik . (95) 

07TG 



Thus the momentum tensor 



1 

T ik = npiU k + VS ik with V = -— — g k , (96) 

07TG 



is not symmetric and the pressure V is negative. 



The Hamiltonian ( 70 I (generator of time translations), expressed in terms of the momentum 

J d 3 xT tt {x,t) (97) 



density, takes the form 

H 



with 

1 2 

T tt = nmpi - —^g^. (98) 

87TG 

The energy density T u , of course, obeys the local conservation law 

d t Tu = -d k T tk (99) 

with 

T tk = u k (T tt + V) . 

If we identify, as usual, the space-time component of the energy-momentum tensor with 
the momentum density 

T a = npi (100) 

we may conclude, that the energy-momentum tensor T a ,p with a, /3 G (t, k) obeys the standard 
conservation law 

dtT 0li + d k T^ k = 0. (101) 
But, because T a @ is not symmetric, we note that 

d t T t p + 8 k T k0 £ 0. (102) 

We remark that such a nonsymmetric energy-momentum tensor was introduced for the 
first time in 1947 by Weyssenhoff and Raabe [23] . 

5.2 Zero mass Galilean symmetry 

Galilean boosts are generated by 

Ki = tPi - ^ j d 3 x 9i , (103) 
which is a conserved quantity because the EOM for the gravitational field gt is given, due to 



I73J and |85|, by 

gt = npi - d k (u k gi). (104) 
The PB of the boost generator Ki and the momentum Pj vanishes 

[Ki, Pj] = (105) 
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as is easily see from the fact that Pi is the generator of space translations 
[Pi, A(x,t)] = di A(x,t) for any A. 



(106) 



The vanishing of the right hand side of ( 105 1 tells us that the total mass of the darkon fluid 
is zero. 



5.3 Angular momentum 



The continuum generalization of expression ( 19 1 for the one-particle angular momentum is 
given by 

Li = e lkl n J d 3 £(x k (£t)pi(£t) + x k (Ct)qi{i,t)). (107) 
By using the identity ( |41[ | and by passing to the Eulerian variables we obtain 

Li = dikl 



J il'.r \ xj,npi - 737=; ''/ !H 



4ttG 



(108) 



5.4 Anisotropic scaling 



Anisotropic scaling, characterized by a dynamical exponent z / 1, implies that the time and 
space variables transform differently w.r.t. dilations 



t — > t* = At, x — > x* = A z x, with A > 0. 



(109) 



In section 2.1 we have shown that the dynamics of free massless particles allows for an 
implementation of dilations with an arbitrary z. To see what happens in the case of a darkon 



fluid we consider a continuum generalization of ( 22 1 of the generator of dilations D 



D = tH — — I d 3 xnpiXi 



AtyG 



d xmgt, 



(110) 



where in the last term we have used the relation (41 1 between the field q and the gravitational 



field g. Taking now the time derivative of D QllOp and by using the EOMs (|73|) we obtain 

5 



dt SttG 



-3 + 



,3 2 

d xg t . 



(Ill) 



Thus D is conserved if z takes the value z = |. It may easily be checked that with z = | 



all PBs between D and the elements of Go are in agreement with the PB-relations (21 1 



The dilations are implemented by the following transformation of the fluid fields A 

A(x,t) -> A*(x,t) = A(x*,t*), (112) 
where za denotes the scale dimension of A. . For an infinitesimal dilation we get 

[A(x,t), D] = (tdt + ^(x k d k + z A ^j A(x,t), (113) 
which leads, due to the explicit form of D ( 110 1 and the PBs (761, ( 86)88 1, to the following 



values for the za- 



Z n = 3, z = 1, 



(114) 



and therefore, due to ( |85[ ), z Pi =1. 

Let us finally prove that our EOMs are not invariant with respect to expansions. We prove 
this by demonstrating the contradiction. So consider the expansions: 

1 



f 



t 



n(t)t, 



n(t)*x with Q(t) 



Kt 



(115) 
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Suppose that we would be able to implement such expansions. Then there would exist a 
generator K of expansions which obeys, in combination with H and D an 0(2, 1) algebra. 
On the other hand, as is well known [18], the PB-algebra between the elements of Go and the 
0(2, 1) algebra would not close for 2 7^ 2. The PB between K and the generator of boosts Ki 
would then imply the existence of a conserved generator of accelerations F, 



[K, Ki] = 




(116) 



But we know from the discussion given in section 3 that the relation ( |41[ ) between the 
field q and the gravitational field g destroys the invariance of our EOMs w.r.t. arbitrary time 
dependent translations (gauge transformations), in particular w.r.t. to constant accelerations. 
This contradicts the initial supposition of the existence of expansions as a symmetry of the 
EOMs. 

Note that the coupling of our darkon fluid to baryonic matter (see 3.4) will break dilation 
symmetry. 



6 Darkon fluid cosmology 

We consider the darkon fluid as a model for the dark sector of the Universe. Hence we look 
for solutions of the fluid EOMs satisfying the cosmological principle (the Universe is supposed 
to be isotropic and homogeneous on large scales). The relevant cosmological EOMs are listed 
in subsection 6.1. In subsection 6.2 we derive and discuss the corresponding Hamiltonian 
dynamics. We start with the three cosmological EOMs in lst-order form, derive the two 
constants of motion and construct out of them that Hamiltonian which has the correct scaling 
dimension. Unique expressions for the corresponding PBs are obtained if we consider explicitly 
the implementation of dilation symmetry. For that purpose it is convenient to introduce new 
independent variables with either scale dimension zero or one. The PB-algebra possesses one 
Casimir which can be used to obtain a foliation of the three-dimensional phase space into an 
infinite number of two-dimensional phase spaces. Within the two-dimensional phase space we 
get a completely integrable system. We are able to construct a corresponding Lagrangian. 

Due to the appearance of a constraint it is not possible to quantize the arising Hamiltonian 
by applying the Dirac quantization rule. An alternative way seems to be the afflne quantization 
(see [26]). 

In subsection 6.3 we solve the cosmological EOMs with the cosmological redshift z as the 
independent variable. By choosing some integration constants to be positive our model shows 
qualitatively the transition from a deceleration to an acceleration phase of the Universe. Fixing 
two integration constants by observational values we obtain a prediction for the cosmological 
history in general agreement with the existing data of the Hubble parameter. 



6.1 Cosmological EOMs 

The realization of the cosmological principle within the Lagrangian formulation leads to the 
following structure for the field variables 

Xi&t) = a(t)$i, (117) 

n(x,t) = n(t), (118) 
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where a(t) is the cosmological scale factor, and 

n(t) 



gdx,t) = g(a(t))xi, n = n(0). 

n 



(119) 



The split of the r.h.s. of (1191 into two factors and defining 5 as a function of a(t) are a 
matter of convenience. If we next put all this into the original EOMs we obtain the following 
cosmological EOMs 



• From the definition of n (see ( 36 1 ) we get 



i(t) = 



From the EOM (42 1 and using ( 120 \ we get 



o(t) = 



n 



g( a ) 

a*(t)' 



(120) 



(121) 



which is a Friedmann-like equation. 



• From the EOM for gt (43 1 and the definition of 0< (see ( 47 () we get 



9( a (t)) = ^T77V> K ± = const., 



where K\ is related to #i by 



m = -— 6 

n 



(122) 



(123) 



We note that possible derivative terms in the gravitational coupling ( 29 \ give no contribution 



to the cosmological EOMs. Taking into account, additionally, baryonic matter (see 3.4) will 
not change the cosmological EOMs. In such a case g(a) would be changed by an additive 
constant (see [T5]V 



6.2 Hamiltonian dynamics 

The cosmological EOMs \l2l\ and (Il22b , when written in l st -order form, take the form 



a = b, 6=4, g = ^. (124) 

So we have to consider a dynamical system in a three-dimensional space. This finite 
dimensional space, called usually mini-superspace, will be the phase space for our Hamiltonian 



dynamics. The dynamical system (1241 possesses two constants of motion 

2 



Q 2 = bK x 



and 



D a 



(125) 



(126) 



In order to construct Hamiltonian dynamics for the dynamical system (124 1 we have to 



find a Hamiltonian H and the Poisson-brackets (PBs) [x a , xp\ for x a £ (a, b, g) such, that the 
EOMs (11241 become the Hamiltonian EOMs 



H] = U, (127) 

where the PBs are required to fulfill the Jacobi-identities and H has to be a function of the 
two constants of motion Q2 and Q3 only. To reduce the arbitrariness in the construction of 
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H we note that the EOMs (124 1 are dilation invariant, i.e. if x a are solutions of (124 1 then 
by the transformation 



x* a (t) = X Za x a (\i), A > 



(128) 



with z a = — §, 2i = |, z g = | we obtain another set of solutions. For the scaling exponents 
of the constants of motion we obtain 



2 3 

zq 2 = -, and zq 3 = -. 



(129) 



Let D be the conserved generator of dilations. Then we get from ( 127 1 and ( 1281 that 



[^d) — fat ~t~ Z a X a — Tja 

(with no summation over a in the middle term !) leading to 

D = Ht- f, with / = H, 

which is equivalent to 



(130) 



(131) 



[H, D] = H. (132) 

However, this expression tells us that we have zh = 1. Thus neither Q2 nor Q3 but Q2Q3 is 
a possible candidate for the Hamiltonian H. So we choose 



H = 



Kf 



(133) 



where the factor K 1 has been inserted in order to get a simple form for H which will be 
obtained when we rewrite it in terms of some new variables which will be introduced below. 



With this H the PBs for the x a have, according to ( 1271, to be chosen such, that 



r 1 dH - * 

[x a , x fi \ - J a 



(134) 



holds. But (1341 is equivalent to 



where 



Vi) 



-ep^s [x-y,xs\ 



and At 



dH 

dx~, 



(135) 
(136) 



The matrix A a p is not invertible (it is anti-symmetric and of odd order) . Therefore ( 135 1 
has no unique solution for the variables yp . Even if we require in addition, the validity of the 
Jacobi identities, the solutions for the yp will not be unique. To render the PBs unique we 



must consider also the implementation of the dilation symmetry. First of all ( 130 1, ( 131 1 may 
be rewritten as 



[/, x a ] 



(137) 



and 



[/, D] = Ht. (138) 
From \138\ we conclude that Zf = . Thus / is a function of the dilation invariant variables 



A 7 ? 



and y = only. Then from z x 



= [/, x] = [y, x 



and ( 137|l we obtain 
df 



dy 
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and 



= [/, y] = [x, y] 



df 

dx ' 



which shows that the PB between x and y has to vanish 

[x, y] = 0. (139) 

2 

As the third variable (of our phase space) besides x and y we choose w = which has 



1. Then, due to Zh = 1 the Hamiltonian H (1331 has the structure 



H = wh(x, y) 



with 



H x ,y) = [v- 



x\y 



(140) 
(141) 



In order to determine the remaining PBs [a;, w] and [y, w] we next rewrite the EOMs ( 124 1 
in terms of the new variables x, y and w 



x — wixy + 3), y = 



2w 



V - 



and w — w 
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(142) 



(1401 



But ( 142 \ are required to be EOMs derivable from the Hamiltonian H. So we get with 

y-D- (143) 



[x,H] = [x,w]h = w(xy + 3) and [y, H] = [y,w]h 



2w 



Solving ( 143 1 for the two PBs gives us 



[x, w] 



w(xy + 3) 
h : 



[y, w] 



xh 



It is easy to check that the PBs (1391 and (1441 satisfy all Jacobi-identities. 
Expressing now the original variables a, b and g in terms of x, y and w 

I _ 3 I 2 | 1 

a — (xw) 5 x, b = Kl (xw) 5 y, g — (xw) 5 

we obtain their PBs 
[a, b] 

and 



H -i ('2 fe2+ 3ff 
5 5 a 



[b,g] = 



[a,g] = 

2 ^ 



_! ^6g + 3iG 

5 5 a 



(144) 



(145) 



(146) 



5 a 2 Q 3 



The equations ( 146 \ can be summarized by one equation 

[x a , Xp] = H'^^faZpXp — fpZ a X a ). 



(147) 



(no summation over a or ft.) The PB-algebra (1471 is identical with a construction given 
by Hojman |24] who starts also with some EOMs, a Hamiltonian and a symmetry. But in 



contrast to |24| our derivation of ( 147 1 shows for the particular case of dilation symmetry the 
uniqueness of the PB-algebra. 

6.2.1 One exact solution 



The EOMs \142\ have a simple power law solution of the form 

1 1 5 



x{t) 



-6, y(t) 



(i 



(t-t ) 



(148) 



which for to = are dilation invariant. We will later show that ( 148 1 governs the behaviour of 



the cosmological solutions at early times. Furthermore the solution ( 148 1 determines a relative 



minimum of the Hamiltonian ( 140 1 at H — 
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6.2.2 Reduced phase space 

In the following we will solve the EOMs < |142| l for x as a function of y. So we reduce the 
dimension of phase space from three to two. To do that we can either use standard analytical 
methods or algebraic methods. The latter shows more clearly the physical nature of this 
reduction. For this reason we will present both methods. 
Analytical method 



From the first two EOMs in ( 142 1 we obtain immediately an ODE of Bernoulli-type for 



x(y) 



where x (y) 



x (y) 



x(xy + 3) 



dx(y) 
dy 



l-2y 

As usual we use the transformation 



?(y) 



to obtain a linear ODE for u{y) 



u (y) + 



3tt(y) 

l-2y 



*(v) 



+ 



1 



2y 



(149) 



(150) 



(151) 



The point y = | is a singular point of the DE ( 151 1. Therefore we obtain one solution u_ (y) 



for y < | and another one u+(y) for y > |. 



w-(y) 



y + k- 



V 



■w+(y) 



y + k+\y 



(152) 



1 

3 •' " \2 

with real-valued constants k+ and Note that the first term in (jl52h represents a particular 



solution of the inhomogeneous DE (1511 whereas the second one represents the most general 



solution of the homogeneous part. The two solutions u± are not connected to each other by 
an analytical continuation in y. This implies that the reduced phase space consists of two 
dynamically disconnected parts with either y < | or with y > h ■ We will show later that the 
validity of the condition 

y > \ (153) 

is necessary in order to explain the accelerating expansion of the Universe. So we will use 
y > \ exclusively in the following. 
Algebraic method 



The PB-algebra (139 1, (144 1 possesses a Casimir 

C = Q3Q2* 



(154) 



where we have used the equivalence between ( 1531 and the condition 



Q2 > 0. (155) 
To prove the Casimir-property of C we observe first that zc = 0. So C can depend only on 



x and y. By using the definition of the Qi, i = 2,3 and the transformation formula (1451 we 
obtain 



C=(y-b 



+ y 



(156) 
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Therefore, due to (139 1, the PBs [x,C] and [y,C] both vanish. To show that [w, C] = 

:ring the chain of relat 

= [H, C] = [w, C] 



too we use ( 139 I again by considering the chain of relations 

dH 



dw 



(157) 



But as |f / we see that [w, C] must vanish. 

This reasoning shows that any function of the independent variables which is a constant of 
motion and has a vanishing scaling exponent is a Casimir of our PB-algebra. But there exists 



no other function, functionally independent from (1541, which possesses this property. Now 
we use the fact that the equation 

C(x,y) = k, ke R 1 , (158) 
defines, for different values of k, a foliation of the three-dimensional phase space (cp [25]). So 



we may solve ( 154 1 for x at fixed k and obtain 

x(k V) = 1 

\-V + k[y-\) 

in accordance with ( 152 1 if we identify k — k+. 



V > 



1 



(159) 



This way we see that our system is described by a Hamiltonian dynamics on a two- 
dimensional phase space with variables w and y and the Hamiltonian 



H = w ( y — 



l + x(k,y) y 



2 J \ ' V 3 
where for the last equality we have used |159j), and the PB 

r i 2w 

[w, y\ — ; 

kx 2 (k,y) (y- 1)5 

leading together to the following Hamiltonian EOM 



= kwx(k,y) ( y- 



1 \ 2 



y = 



2w 
c(k,y) 



x(k,y) 



(160) 



(161) 



(162) 



The dynamical system ( 162 \ is completely integrable. To show this we make use of the 



existence of a second constant of motion besides the Casimir (see (1251, ( 126 ) ) . In this proof 



we may use the Hamiltonian H ( 160 1. Let us put H — E, i5=const. Then we can solve ( 160 \ 
for w 



w(E,k,y) = 



E 



kx(k,y) (y- 1)1 



(163) 



Inserting ( 163 1 into the first eq. in ( 162 1 we finally obtain y(t) by a quadrature. 
Generator of dilations 

We can also give an explicit expression for the generator of dilations in two-dimensional 

'By 



phase space. From / = H ( 131 1 we get y^- = H and so 



F(y)- 



(164) 



Lagrangian 



The EOMs ( 162 1 are Euler-Lagrange EOM derived from the Lagrangian 



yF(y) log™ - H(w,y), 



(165) 
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where H(w,y) is given by (1601. To prove this it is sufficient to note that 



U(y,w) = F(y) \ogw 



(166) 



is the momentum canonically conjugate to y. Indeed, by using ( 161 1 we obtain 

[y, Tl(y,w)} = 1. (167) 



We can now use ( 166 \ to eliminate w from the Hamiltonian (j 160p to obtain 

1\ ' 



H = kx(k,y) y 



cxp 



{nv)) 



(168) 



With |167j and ([168]) we have finally Hamiltonian dynamics in a two-dimensional phase space 
expressed in terms of canonical variables y and II. 
How to quantize this system? 



We would like to quantize the system (167 1 and (168 1 with the constraint (1531. This 



cannot be done if we try to follow Dirac's rule of going from the classical PB (167 1 to the 
quantum commutator 

[y, ft] = ih. (169) 



The reason for this is that, for self-adjoint operators y and II (1691 has, up to unitary 
equivalence, only one irreducible representation given by 

f (— oo < y < oo) (cp |26|). So the spectrum of y covers the whole real axis. To 
overcome this difficulty we follow J. Klauder (see [5B]) and consider instead of the canonical 



quantization ( 169 1 the affine quantization given by 

[q, £>] = ihq 



with q = y — \ and 



D EE Uilq + qfl) 



(170) 



and choose the irreducible representation for q > 0. Details of such a quantization procedure 
still have to be worked out. 



6.3 Solution of the EOMs 

Let us first put the two constants of motion Qi,z equal to numbers 7^2,3- Next we solve the 
defining equation of K2 ( |125[ l for a — b and obtain 

a = ± U + £) . (171) 



So, with K2 > 0, equivalent to condition (1531, we must have 

Ki > (172) 
in order to describe an expanding Universe. 



Using the definition of K3 ( 126 1 we obtain (we consider now g as a function of the redshift 

2, 



(1 + z) 1 se a), a cubic equation for g(z): 



" 3 '1 + 2 



T+^M 1 -!^) (173) 
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where the transition redshift z t is defined by 

1 



(174) 



1 + zt K s 

Let us look now at (|173[) for K2 > and K3 > together with the Friedmann-like equation 



(see 1 121 1) 



g{z){l + zf 



(175) 



For z > z t we have g < and so, due to ( 175 1 a < 0. So, for a > zt we are in the 
deceleration phase of the early Universe. On the other hand, for z < z t we have g > and 



then, due to ( 175 1, a > 0. So, for z < zt we are in the acceleration phase of the late Universe. 
To get these results we have to require that iG,3 > 0. 
To see this note that 

• From (11741 we conclude that K s > 0. 



• Let z t be a zero of g. Then differentiating ( 173 1 at z = z% shows that ^g (zt) < 0. But 



the transition from the decelerating to the accelerating phase requires g (z t ) < 0. . Thus 



we must have & > 0. 

«3 



Next, to simplify ( 173 1 we introduce 

g(z) = K-^g(z). 



Then ( 173 1 becomes (we recall that k — KsK 2 

f + S = k[ J 



1+2 



1 + Zt 

where g now depends on the two parameters z t and k only. 



(176) 



(177) 



Using ( 171 1 the Hubble parameter H[z) = -(z) may be expressed in terms of g 



H(z) 



K2 

Ki 



1 + 



(1 + z). 



(178) 



However, it is more convenient to consider instead of H (z) the normalized Hubble param- 



eter 



for which we find 



h(z) = 



H(z) 
Ho ' 



h(z) 



1 + 



(179) 



(180) 



So, if we know Ho, the cosmological history is determined in our model completely by the 
two parameters z t and k. 

To obtain k in terms of observable quantities we consider the deceleration parameter 



a h (z) . 



q(z) = -— = — -1| 1 + ;) - I. 
which, when expressed in terms of g, Zt and k takes the form: 



q{z) = - k ^^g{z)[l + 

1 T Zt 



(181) 



(182) 



Thus the value of k can be determined from the knowledge of the present deceleration 
parameter go — q(0) and of zt- 
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Remark: The expressions for g(z), h(z) and q(z) in our paper Q2] can be obtained from the 
present ones by the substitutions g(z) — > — v6<?(z) and k — > yf&k. Furthermore, the constants 
K\ and no are related to the constants D and /3 from 15 by the relations K\ = —3D/3 and 



3D 

4-7T ' 



Asymptotic results for z>>zi 



From (177 1 we obtain 



m ~ - (^rr^ ) • (183) 

Thus we obtain a power law behaviour for 

ft,(z) oc zs, (184) 

which shows that asymptotically our cosmological solutions are scale invariant. On the other 
hand the deceleration parameter reaches asymptotically a constant value 

q(z) ~ |. (185) 



The asymptotic behaviour ( 184 1 differs from the predictions of the ACDM-model 



hACDAt(z) ~ 25 (i 8 6) 



for the matter dominated epoch. But we expect a change of our prediction ( 184 1 if we generalise 
our model to the one that exhibits general covariance (cp. Section 9). 

6.4 Distances 

We must decide how to relate observational data of distances (e.g. of supernovae) to the 
Hubble parameter H(z). In standard FLRW cosmology we obtain, in terms of comoving 
coordinates for a flat Universe and for a radial light ray (see section 1.2 in [57]) the coordinate 
distance 

« w - C I w-y (187) 

In any relativistic theory of gravitation which possesses a metric the limit of an isotropic, 
homogeneous and flat Universe is described by the FLRW metric given in comoving coordinates 
by 

ds 2 = a{t)(d£) 2 - c(d£) 2 . (188) 



Then eq. ( 187 1 arises from the requirement that light propagates along a null geodesic 

ds = 0. (189) 

Note that our darkon fluid model is neither a relativistic nor a metric theory. However, it may 
be understood as a limiting case of such a theory as shown in section 9. On the other hand, 
in a nonrelativistic theory we do not have an expanding space but only moving particles in a 
static space. But both points of view lead to the same observable effects as has been already 
noted by Milne in 1934 [19] . We need only to assume l that, relative to a particular observer O 
the speed of light is a constant c independent of the light source and that the classical formula 
for the Doppler effect is valid 1 (McCrea, page 351 in [28]; see also section 5.2 in [29] ). So both 



points of view allow us to use the expression ( 1871 for the coordinate distance in our model 



2G 



The dimensionless coordinate distance £o, denned by 



= ^tiz) = r ^ (wo) 

c J h(z ) 



may easily be computed in terms ol the two parameters Zt and k by using (180 1 with (1771 
leading to (cp. Appendix A in |15| ) 



&(*) = 6 (l + ^M) £ a t log(x - xOlfgJ. (191) 



%2 A 

2 

where 

a t = {{xi - Xi+i)(xi - Xi-^y 1 (192) 
i = 1, 2, 3 and cyclic permutations and the a;j are the roots of the cubic equation: 

a: 3 + 6 (as - ft) = (193) 

given by 

x\ = u+ + w_, a; 2 = 1 iv 3, 13 = ^2 

with 

v± = (3fc ± (8+ (3fc) 2 )5)i 
In terms of ^0 we obtain, as usual, the dimensionless luminosity distance Dl(z) 

D L (z) = (l + ^ (z), (194) 

the related distance modulus fi(z) 

fi(z) = 51og 10 D L (z) - 51ogfl"o|ioofcm S -iM P c-i + 42 -38 (195) 
and the angular diameter distance Da{z) 

D A (z) = D L (z){l + z)- 2 . (196) 

6.5 Predictions versus observations 

In subsection 6.3 we have shown already that our model explains qualitatively the transition 
from a deceleration phase of the early Universe to an acceleration phase of the late Universe 
if we choose the constants Ki[i = 1, 2, 3) to be positive. In this subsection we will show that 

• the quantitative predictions for the Hubble parameter H(z) are in a general agreement 
with the observational data, 

• the distance moduli for GRBs seem to show a peak at large z for the deceleration pa- 
rameter q(z) in agreement with our predictions, 

• our unified model for the dark sector is unable to predict the present dark matter con- 
tribution QmO and, therefore, it can neither explain the BAO data nor the data on the 
CMB-shift parameter by using the present models for them. 

At present we do not know what the changes for the prediction of the Hubble parameter 



( 180 1 would look like had we generalized our model to a general covariant theory presented in 



section 9. Therefore we have not performed, so far, any 'least squares fits to the data. 
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-1 -0.5 



Figure 1: Comparison of experimental data from [32], |33| . |34j for H(z) measured in units of km 
s _1 Mpc -1 with our predictions for zt — 0.71, qo = —0.57 and Hq = 70.5 (dotted line) and a 
ACDM fit for H = 73.8 and O roi0 = 0.27 taken from [33] (dashed line). 



6.5.1 Hubble parameter 

In fig 1 we compare the best out of three predictions of our model for h(z) given in [15] with 
observed Hubble parameters which we adopt from table 1 of [30]. In addition we show the 
ACDM-fit to these data taken from [38] 

In [15] we used z t = 0.71 and qo = —0.57 taken from a ACDM- fit to some data given 
in [31J . Let us point out that these values for zt resp. qo are in agreement with very recent 
model independent constraints on cosmological parameters [5S], [3B], [37]. For our prediction 
we have fixed the normalization in fig. 1 by choosing Ho = 70.5km s -1 Mpc -1 in accordance 
with the value used very recently in |36| . The corresponding values for the three constants Ki 
turn out to be 

Ki = 2.361 x 10 5 (km s' 1 Mpc' 1 ) 3 

K 2 = 1.263 x W^kms^Mpc- 1 ) 4 
K s = 9.530 x 10 10 (km s' 1 Mpc' 1 ) 6 
leading to k — 2.123. 



= 8.033 x 10" 5 s 3 

= 1.392 x 10" 71 s" 4 
= 1.103 x 10~ 106 s - ' 



6.5.2 Remarks on distance moduli data 

The most recent supernovae (SN la) distance moduli data are those of the Union 2.1 compi- 
lation 39 . But we do not see any necessity to give a detailed comparison of these data with 
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the predictions of our model (see eq. (195)) because it has been claimed very recently by O. 
Farooq et al. [3D] that the recent 21 Hubble parameter versus redshift data points (see fig. 
1) give constraints on DE models which are almost as restrictive as those from the Union 2.1 
compilation. In addition to the SN la data we have for larger z values up to z — 8.1 distance 
moduli data from Gamma- Ray Bursts (GRBs) at hand (see table II in [ID])- Using such 
data one obtains the tentative result that the deceleration parameter q(z) shows a maximum 
roughly at z ~ 1.6 (41] in qualitative agreement with the predictions from our model (see fig. 
2 in [15]). 

6.5.3 Remarks on BAOs and the CMB-shift parameter 

Baryon acoustic oscillations in the early Universe are imprinted as peaks in the late time 
matter power spectrum (cp. 42 ]). But the physical interpretation, mainly of the first peak 
position, depends heavily on the present baryonic and dark matter contents Qto resp. f2 m o of 
the Universe (see fig. 4 in [12] ■ This issue is usually ignored in the data analyses performed 
in the literature 43 . But our model, as a unified model for the dark sector of the Universe, is 
unable to predict f2 m o (for more details cp. Section 6.6) and therefore we cannot at present 
analyse these peaks. The same argument holds for the CMB-shift parameter (cp. 44 ). This 
issue has been discussed already in our recent paper [15] . 

6.6 Comparison with alternative models for dark energy 

In the following we will sketch how to perform a comparison of our model with alternative 
models of the dark energy. Any details will be left for the future. 

The Friedmann equations, derived from GR including a cosmological 'constant' A, are for 
a flat Universe (see e.g. [27]) 

u2 _ 8nGp + A 



3 3 

a 4-kG A , „ . 

a =-^-0> + 3p)+ 3. (198) 

where p is the total energy density of the Universe (sum of the matter (M) and dynamical 

dark energy {DE) if we neglect radiation) and, if we assume that matter is described by 

pressure less dust, p is the (supposed negative) pressure of DE 

P = PM + PDE, P = PDE < 0. (199) 

Let us now assume that the first Friedmann equation 11971 holds. Then the second equation 



( 198 1 can be replaced, as is well known, by the continuity equation 



d_ f _K 
dt \8n( 

Note that if we allow a weak time-dependence of Newton's gravitational constant G the 



dtK87vG +P)+^(p + P )= 0. (200) 



last equation (2001 would have to be generalised (see [45], [46]). However, this case will not 
be discussed here, instead we refer the reader to |46] and the literature cited therein. 

At this stage it is important to contrast here two possibilities for the description of the 
'dark energy' [45] : 
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Some kind of fluid described by an equation of state (EoS) W(z) 



W(z) 



P 

pDE 



• A (possibly) time-dependent A with a fixed EoS 

A 



PA 



= -PA- 



Now our Friedmann-like equations are (see (175) and |178l) 



(201) 



(202) 



Kj g(z) (1 + zf 



and 



I<2 

Ki 



1 + 



(i + -) 2 



(203) 



(204) 



Comparing the corresponding right hand sides with ( |197|198j ) we obtain two equations for 
the effective functions p, pa and p 



P- PA 



P + PA 
1 

~8^G 



3H 2 

2 



3 

8vrG 



1 + 



I<2 
2 



1 + 



9\A 



1 



(1+*) 



■Mi g( Z )(l + zf 



(205) 



(206) 



2 ) v ' 4ttG 
In the first scenario we have to put pa — and arrive at two equations for the three 
functions pm, Pde and p. If we suppose that we have no interaction between CDM and DE 
we have, in addition, local conservation of pM 



with the solution 



PM + 3Hp M = 



PM = pM, o(l + z) 3 



(207) 



(208) 



From ( 205 I, ( 206 1 and ( 208 1 we easily get an equation for the effective EoS which, however, 



contains as a free parameter the present value of the energy density of matter (CDM + 
baryonic matter), pu.o, which is not determined in our unified model for the dark sector of 
the Universe, i.e. pm,o is a free parameter in our model (for a very detailed discussion of 
this issue see [47]). Therefore our model, as any other unified dark fluid model, is not able to 
predict an EoS for DE. However, on the other hand, we can define an EoS Wd for the darkon 
fluid as a whole 

2Kig(z)(l + z) 

P 



2 



(209) 



3(f) ^1 + 

Whether Wd for some z-region shows a phantom- like behaviour (wn < — 1) will be left for 
future considerations. 

In the second scenario we have to put pDE — p = which corresponds to model II in |46j . 
We decompose pu into its DM and B parts 



PM = PDM + PB 



(210) 



and take into account the local conservation of ps which leads, in the cosmological limit, to 
(cp. p08t) 



pB = PS,o(l + z) 



(211) 



Note that Pb,o, in contradistinction to pm,o, is directly observable. Then (205 1, (2061 lead 



to predictions for pdm and pA for a given set of constants Ki (i — 1,2,3). More details and 
their comparison with the results communicated in [46] will be left for the future. 
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7 Modeling dark matter halos by a steady state darkon 
fluid 

We start with the darkon fluid EOMs in the Eulerian formulation as given in section 3. By 
taking their time-independent form (steady state equations) we will derive a nonlinear ordinary 
differential equation (ODE) for the gravitational potential in the spherically symmetric case. 
The solutions of this ODE may serve as a model for the dark matter (DM) halos of galaxies. In 
particular we consider DM-dominated galaxies in which case the halo dominates the galactic 
gravitational potential and therefore determines the rotational velocity of stars (rotation curves 
(RCs)). 

In subsection 7.1 we derive the nonlinear ODE for the gravitational potential in the spher- 
ically symmetric case and the behaviour of their solutions for large distances from the centre 
of the halo. Then, in subsection 7.2, we present the results of the numerical solutions of this 
ODE and compare them with the main characteristics of the observational data for the RCs 
of DM-dominated galaxies. 

7.1 Equation for the gravitational potential 

From the darkon fluid EOMs ( |73[ ) we obtain immediately the following steady state equations: 
Steady state equation for the gravitational field 

u k d k (H) = 9 —d lUk - 6i, (212) 
V n / n 

where the auxiliary field 8i is given by the solution of 

UkdkOi + 6 k diU k = 0. (213) 

Euler equation: 

UkdkUi = Qi (214) 

Continuity equation: 

d k {u k n) = 0. (215) 

Next we consider the halos as spherically symmetric objects. Then the vector fields ap- 
pearing m ( |212|215[ ) take the following form: 



u k (x) = —u(r), g k (x) = - — <j> (r), (216) 
r r 

8 k (S) = 5»0(r), 
r 



where cf> is the gravitational potential. By using 1 216 1 in the steady state equations we obtain 



from (2121 



an 



« — =0, (217) 



where 0, obtained by the integration of ( |213[ ), is given by 

^ _ 1 7 



(218) 



2 au{r) 

For convenience we have chosen to write the integration constant as — 3- . Euler's eq. (214 1 
leads to the Bernoulli eq. 

\u = -<t> (219) 
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and the solution of the continuity eq. gives 



a 



(220) 



where a is an integration constant. 

As in the Lagrangian formulation the vector field 8 is a conserved quantity we may wonder 
whether it is possible to connect the integration constants a and 7 appearing here and K\ 



appearing in the cosmological context (see ( 123 ) ) . 
We find, in our case, 



9,(1) = (221) 



This expression has a different £ dependence from that of the r.h.s. of (123 1. So, unfortu- 
nately, there is no relation between a and 7 on one hand and K\ on the other. 

We note that the solution (12201) produces a singularity at the origin on the r.h.s. of the 



continuity eq. (2151. This singularity might be interpreted as the central black hole of the 
galaxy (see [48] for the existence of a central black hole in LSB-galaxies). But a black hole 
describes a sink for any kind of matter. So we get a < and therefore 

u < 0. (222) 



Inserting all these expressions into (2171 we obtain a ODE 

^V)' = ^(-20)- f (223) 

which may be understood as the Poisson eq. for a self consistently determined effective DM 
mass density p 

4nGp(r) = JL(-2*)-i, (224) 
Next we switch over to a dimensionless potential <j> by the substitution 

4> -» ji<f>. (225) 

From the foregoing equations we easily conclude that 7 has the dimension m 5 s -5 as it should 
be. 

So, altogether, we end up with a nonlinear, non-autonomous ODE for <f> 

(rV)' = |(-2^)-i, (226) 
in which <j> is nonpositive and satisfies the natural boundary conditions 

4>(oo) = and 0(1) = -|, 77 > 0. (227) 



We note that the ODE (2261 and the first boundary condition (<p(oo) = 0) are invariant 
w.r.t. to the scale transformation 

4>(r) -> 4>(r) = cf>(\r), A > 0. (228) 

Thus the second boundary condition is necessary to make the solution unique. 

To simplify ( |226[ ) we introduce, instead of <j> the function ip — —2rcj> which then has to 
satisfy the ODE 

ip" = -r-5 ip~Hr), (229) 
whose positive solutions have to be determined. 
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Remark: Eq. ( 229 1 resembles the Thomas- Fermi equation which arises in the self- 



consistent determination of the potential within a neutral atom. But there are two important 



differences between these two problems. Besides the negative sign in ( 229 1 , the function on 



its r.h.s. is reciprocal to the one appearing in the Thomas-Fermi equation. 

To determine the required solutions of ( |229| l let us consider first ( |229[ ) with the boundary 
conditions 

= 77, 77 > and tp' = /3. (230) 

In this case we can use the following lemma, which is a particular case of a more general 
theorem derived by S. Taliaferro (see theorem 2.1 in I49|) 
Lemma 



Let ip(r) be a positive solution of (229 1 satisfying (230 1 for fixed rj > 
and an arbitrary real number (3. Then there exists r%, satisfying 1 < T2 < 00 such that 
lim r ^ r2 tp(r) — and ip(r) cannot be continued past r2 as a twice continuously differentiable 



solution of ( 229 1 



To use this lemma to the study of properties of required solutions of ( 226 \ it is advantageous 
to perform first the transformation 



x = logr, ip(r) — > <p(x) 



-X 1 1 X 



(231) 



Then ( 229 1 with the boundary conditions ( 230 1 gets transformed into the following au- 

(232) 



tonomous ODE 



p (x) + ip (x) = — ip 2 (x) 
with the boundary conditions 

((3(0) — a, a > and tp (0) 

where a and b are dimensionless numbers. 



b,beB} 



(233) 



By (231 1 the lemma may be translated to the following statement about ip: For fixed a 

<p(x) = 0, where X2 = logr2, and tp cannot be continued past 



and arbitrary b we have lim^ 



X2- Eqs. (232 I and (2331 are together equivalent to the integral equation 

1(8). 



(p(x) 

Using the relation 



6 1 



c x 



ds (1 



9 



^V(logr) 



(234) 



(235) 



( 234 1 may be rewritten as an integral equation for 



('■) 



1 

1 

2 



ds(-2(f>(s)Y 



ds 
s 



—a — b 



ds 



(236) 



+ 



ds(-2cf>(s)y 



which is cast in such a form that the first line contains the usual integral of the Poisson 
equation, the second line a constant and the third one a term being proportional to 1 . 



But the latter is in conflict with the original ODE (2261. The reason is that the transfor- 



mation (235 1 becomes singular at r = 0. Therefore the coefficient of i has to vanish which 

-r 

J — c 



fixes the value of b at 
b 



I 

Jo 



ds{-2<t>(s))~ 



dse" ip 2 (s) 



(237) 
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But the lemma given above holds for any b, and so in particular also for b — b c . Thus, 



due to the relation (235 I, it holds also for <f>(r) ie. lim,-_j. r2 <J>(r) = and so <j>(r) cannot be 
continued past ri- 



Next we return to looking at solutions ip of ( 232 1 . We know from the results given above, 
that 

• we have to find a positive solution tp at least for < x < X2 with ip(x2) = 0. 

• we have to vary b until it reaches a value b c which corresponds to 

r° 3 

b c = - / dse <p~2(s) (238) 

J — oo 

leading to 

/ f X 3 

ip < -e~ x / dse s ~ x <p~i(s) < for all b < b c . (239) 



To relate the solutions of (2321 to the RCs v(r) of a DM-dominated galaxy we recall that 
v(r) is determined from the equality of the centripetal acceleration and the gravitational force 
acting on a star in circular motion [50] 

^ = -yh'ir), (240) 



where the factor 75 arises due to the dimensionlessness of cj> after the substitution ( 225 1. Then 



using ( 235 1 we can determine v in terms of the derivative of ip 



v 2 (r) = --jiip (x)\ x= i ogr . (241) 

7.2 Predictions versus observations 

To begin with we would like to discuss how to fix, at least in principle, scales for distances and 



velocities. Due to the dilation symmetry of our model (see ( 228 () we have no intrinsic scale 
for distances. A scale for velocities may be chosen by fixing the value of the dimensionfull 
constant 7. However, 7, being an integration constant, is an extrinsic parameter. So our 
model does not possess any intrinsic scale for distances or for velocities. To get intrinsic scales 
we must introduce dilation symmetry breaking terms into the Lagrangian. We can do this 
either by taking into account the coupling to baryonic matter and/or by adding within the 
gravitational coupling ( |26[ ) higher order derivative terms as indicated in ( |29[ ). 

So, for the present form of our model, we get the required scales by fixing two extrinsic 
parameters. For velocities we have to fix 7. For distances we have to fix the scale parameter 



A in (228 1 or, equivalently, put x — in (233 1 which fixes r = R where R is the chosen unit 



of length. Thus, to fit observational data for the RCs of any given galaxy, we would have to 



perform a best fit for 7, R and the dimensionless 'initial' value a in ( 233 1. However, as we view 
the present model only as a first building block for a new theory (cp. section 9) we think it is 
premature to perform any 'least square fits' to the data. Thus we have only checked whether 
there is a general qualitative agreement of the predictions of our model with the data. We 
have done this by fitting the extrinsic parameters 7 and R to one particular galaxy at a fixed 
value of a. Note that the choice of a fixes 6 C . 

More detailed discussion has been postponed until we have generalised our model to a 
general covariant one as suggested in section 9. 
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Figure 2: Numerical solution of (232) for a = 10 and b = —0.03 




Figure 3: Numerical solution of (232 1 for a = 10 and b = b c (thick line) compared with our 



approximation ( 243 ) (thin line) 



We have solved numerically the nonlinear differential equation (2321 for <p(x) by choosing 



the initial conditions (233 1 with a = 10 and b < —0.03. The main result is the abrupt change 



of the character of the graphs for ip(x) at negative a;-values at the critical 6-value b c given 
approximately by 

b c = -0.031475314466, (242) 

For b > b c the graph of f(x) crosses the a;-axis at some Xo < where xq is an increasing 
function of b. We believe, although this is numerically hard to verify, that xo — > —oo for 
b decreasing to b c . For x < xo the solution becomes complex valued. This behaviour is 
illustrated in fig.2 (6 = -0.03) and fig.3 (6 = b c ). 

On the other hand for b < b c the graph of ip(x) starts to rise rapidly at some x < and 
goes to infinity. This behaviour is clear from fig. 4 (b = —0.032). 

Furthermore, from the plot of cp (x) for b = —0.032 (fig. 5) we see that tp (x) < for all 

x e R 1 . 

The behaviour of the graphs of ip and ip for b > b c resp. b < b c as shown in figs. 2 — 5 



strongly suggests that the b c obtained numerically in (2421 is, modulo numerical fine tuning 



indeed identical to the expression for b c derived in section 7.1 as the correct b- value for modeling 
halos by means of the solution of the boundary value problem ( 232|233 1 for a — 10. 
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Figure 4: Numerical solution of (232 1 for a = 10 and b = -0.032 




Figure 5: Numerically obtained <p for a — 10 and b — —0.032 
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Figure 6: Numerically determined plot of the darkon fluid contribution to the rotational velocity 
v(r) (in fcms -1 ) at small radii (in kpc) for b = b c , 75 = 615fcms _1 and R = 4x 10 n fcpc. Compared 
to the RC of UGC 8490 (data from [56]) 



Hence we have to compare qualitatively the RCs calculated for b = b c with observational 
data. For that we plot in figs. 6-7 the rotational velocities v(r) determined from eq. (241). In 
fig 6. we show the behaviour of v(f) at small radii with the extrinsic parameters 75 , resp. R 
fitted, for fixed a = 10 to the RC of the dwarf galaxy UGC 8490 (data taken from [56]). Fig 
7. shows the global behaviour of v(r) in dimensionless units. 

The graphs for the darkon fluid contribution to the RCs given in figs. 6-7 show the 
following qualitative characteristics when compared to observational data on DM haloes (for 
a very recent review see |51p 

• They start at a radius r\ ~ and rapidly increase reaching v(r) ~ 76fcms _1 at r ~ Akpc. 
(fig. 6). 

• Beyond this increase our model predicts a rather flat, but still very slowly increasing 
behaviour for the RCs until extremely large radii. This is basically in agreement with all 
observations according to which the RCs of a galaxy stay approximately constant after 
attaining a maximum at about 5 kpc [52]. This statement is true, in particular, for the 
DM-dominated galaxies: the LSB-galaxies (cp. the RCs reproduced in [53] and in figs. 
3,4 of [54]) and dwarf galaxies (cp. [55] and the RCs reproduced in [56]). 

The behaviour of <p(x) within this flat-like region can be well approximated by (see fig. 3) 

<p(x) ~ /(x) = (2.5) a4 (c-:r) ' 4 , (243) 
where c ~ 126.5 has been fixed by /(0) = 10. Then we obtain /'(0) ~ -0.03162, which 



is a good approximation for b c (242 1. f(x) is the solution of the ODE 

f{x) = -f-i(x) (244) 



and so is an approximate solution of the ODE ( 232 1 if ip' ~ const 



• They terminate at an extremely large radius roughly given by T2 ~ e 130 with a diverging 
RC as r — > r 2 . 
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Figure 7: Global behaviour of the rotational velocity v{r) for b — b c plotted as a function of 
X = log r 

Of course, this contradicts all observations. But the size of r2 is clearly orders of magni- 
tude larger than the size of the visible part of the Universe: In the FRW model the size 
Ro of the Universe is given by _Ro = 0.5c£o [S3, where to ~ 13.7 Gyr is the age of the 
Universe [50]. Thus we obtain Ro ~ 2.1 Gpc. The choice R = 4.10 14 pc in fig. 6 gives 
r2 ~ 10 71 pc, which is about 62 orders of magnitude larger than Ro. 

We conclude that this highly unphysical behaviour of v{r) as r — > r-z is possibly a conse- 
quence of considering only an isolated halo. We expect an improvement of this situation if we 
embed the halo in the time-dependent darkon fluid (to be done by an appropriate generaliza- 
tion of the considerations in the next section 8). 

8 Influence of the cosmic expansion on binary sys- 
tems 

In the following we consider two astrophysical objects, one of large mass M {e.g. the sun) and 
another one of small mass m <C M {e.g. a planet) embedded in the time dependent darkon 
fluid which, asymptotically behaves like the cosmological solution (see section 6.2). We would 
like to assess the change of Newton's force of attraction between the two bodies due to the 
presence of the darkon fluid. To be specific we approximate this binary system by a point-like 
test particle of mass m which moves in the time-dependent gravitational potential created by 
a point-like mass M and the overall darkon fluid considered in the rest system of the heavy 
mass M. 

Within GR the corresponding problem is to find, as a solution of Einstein's equation, an 
interpolation between the Schwarzschild and the FLRW metric. Then the motion of the test 
particle is the geodesic motion in this metric. A solution to this problem has been given in a 
very recent paper [58] (see also the recent review [59 and the literature quoted therein, which 
covers also the history of this problem starting with the paper by Mc Vittie in 1933 [60]). In 
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the Newtonian limit one obtains for the EOM of the test particle (see |58| ) 

Xi = -Xi ^-Xi. (245) 

a r 1 * 

It turns out that the correction to Newton's law by the present accelerated expansion is 
completely negligible within the solar system (see section III A in [59|). 

In the following we will sketch the ideas how to obtain the correction to Newton's law for 
our model. 

We start with the EOMs given in 3.4 with a gravitational potential 

0(r,t) = Mr,t) ~ (246) 

where cf>D represents the darkon fluid contribution. Then we obtain, in the spherically sym- 
metric case, the EOMs 

w + {uw) = 0, (247) 
u + l -(u 2 ) = (248) 
(r 2 <j> D y + u(r 2 <j>' D )' = w9, (249) 



+ (u6 



0, 



(250) 



where the functions u and 9 are given by (216 1 generalized to the time-dependent case and we 
have defined w = r 2 n. A dot represents dt and a slash d r - ■ Next we have to solve ( 247]250 ) 
with the following boundary conditions: 

• for r — ¥ oo the functions w, u, (f>o and 9 attain their cosmological limits defined in 6.1 



w(r, t) 



r 2 no 



u(r,t) 



,«(*) 
a(t)' 



(251) 



f>o(r,t) = 



r 2 g(a(t)) 



(r,t) = -- 



2 a s (t) ' " V, " J ' n a 2 (t)' 
They are solutions of ( 247|250 \ if a and g satisfy the cosmological equations ( 121 1, ( 122 1. 
• for r — > we require that 

r4> D -t (252) 
leading, if inserted into the EOMs ( 247|250 1 , to the following behavior of u, w, 9 and 
4>d for r — > 

u(r,t) ~ /3r~5 with /3 = ±(2GM)?, (253) 

2 



w(r, i) 
9(r,t) 



7-2 /, 

ri f e [t 



3(3 
2 

— v 



5p 



' 5 (Ufo) t 



2 a 



(254) 
(255) 
(256) 



where the functions f w , fe and the sign of /3 still have to be determined. 

If it turns out that it is not be possible to obtain analytic expressions for these functions one 
should try to find an approximation involving a two-scale perturbation expansion. The small 
spatial scale will be given by the size of the binary system and the large one by the cosmic 
scale. This problem is currently under consideration. 
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9 Outlook for a general covariant theory 

The Poincare algebra treats spatial and time variables on the same footing. Then its enlarge- 
ment by dilations has necessarily a dynamical exponent z = 1 (cp. [61]). But in our model 
we have anisotropic scaling with z = §. Therefore, an enlargement of the Milne gauge group 
( |24[ ) to arbitrary space-time diffeomorphisms (Diff) as the symmetry group for a generalized 
darkon fluid model seems to be forbidden. But the dilation symmetry with z = § was an 
emergent symmetry and not an input of our model. So we should look for a general covariant 
generalization of our model without imposing any form of scale symmetry. 

To do this we have first to restore the Milne gauge invariance for the Eulerian action (|56[) 



by undoing the gauge fixing ( 41 1 leading to 

S = ^ J dtd 3 x (4-wGnq^Dtu - g l ) - 8(d t n + d k (nu k )) + naDtf^j + S flM (257) 
with the spatial metric given by r]ij = Sij and S field is determined by the expression (31). 



Now we propose the following action as the Diff-generalization of (257 1 



Ddiff = J d 4 x^^g (^-nGnq v u x A\u" — 9A\(nu x ) + nau x d\P^j + Seh, (258) 

where we have used the following notation (in units with c = 1): g is the determinant of the 
covariant space-time metric g^ v with the invariant line element 

ds 2 = g^y dx^ dx v (259) 

and A\ denotes the covariant derivative 

Aju" = d x u v + Tl„u°. (260) 

Here V^ a is the corresponding Levi-Civita connection, n is the scalar particle density and u x 
the four-velocity. 6, a and fi are scalar fields. The vector field q M will be deduced, as given 



below, from the Milne-scalar qi . Finally the last term in ( 258 1 , which is the Diff-generalisation 
of the last term in ( |257[ ), is the Einstein-Hilbert action given by 

Seh = -jgj— J d 4 xV=gR, (261) 

where R is the Ricci-scalar. 



To understand (2581 as a Diff-generalization of (257 1 we consider the limiting case 

u°-H, ff^-Hfc*, (ij^ = diag(-l, 1,1,1,)) (262) 

TL except To,, -t -g\ 
which are the leading terms known from the Newtonian weak field limit of General Relativity 
(GR) g^v = ?7 M „ + hpu, \hfiv\ <C 1 an <l an expansion of h^ v in inverse powers of c. Here g 1 is 



the nonrelativistic gravitational field as used in (257 1 



To construct the vector field g M we proceed in two steps: 

• By means of a Lorentz transformation we promote the Milne-scalar qi (i = 1,2,3) to 
Qa (a. — 0,1,2,3), which is a scalar field w.r.t. Diff but a contravariant space-like vector 
field, labeled by a, in tangent space. 

• By means of tetrads Eft (known from teleparallel gravity, cp. 62 ) we transform qa to 
the four-vector field q M 

qa -> <7 M = laEft ■ (263) 
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The tetrads Eft transform like a contravariant vector w.r.t. Diff by changing the Greek 
index and like a four- vector w.r.t. Lorentz transformations by changing the underlined Greek 
index. The tetrads are related to the metric tensor by 

g^ = v<*l E » E »- ( 264 ) 
So the limit t/ M „ — > r)^ corresponds to Eft — > 8ft and so to <J M — s> g M . Therefore we obtain, in 



the limiting case, for the individual terms in (2581 

q v u X A\u" — > qi(D t u — g'), 
A\(nu x ) — > d t n + dk{nu k ), u x d\j3 — > D t f3 



(265) 



in agreement with (2571. 



What about possible scaling properties of the action (258 1? The generators of Diff do not 



form a closed algebra together with the generators of anisotropic scaling transformations 



Xi 



= Xx l 



(266) 



A closed algebra can be obtained only if the rigid transformations ( 266 1 are replaced by 



local ones (Weyl transformations). But the Einstein-Hilbert action (|261[) does not respect the 



Weyl symmetry (for that we would have to replace (261 1 by the Weyl action (cp. [63] and 



the literature cited therein)). So we can only expect, like in Horava gravity [64], to find scale 
invariant solutions in the ultraviolet or in the infrared regime. Any further consequences of 



this new ansatz (2581 for a general covariant darkon fluid have still to be worked out. In 



particular we are interested in the cosmological implications. But due to the larger number of 
gauge fields in the action (instead of the three gravitational fields g % we have now ten fields of 
the metric tensor g^v) we will obtain cosmological EOMs other than the ones studied so far. 
In particular we will obtain a prediction for the Hubble parameter h(z) which would probably 



be different from (1801. Therefore nothing can be said at the moment about the asymptotic 
behaviour of h(z). 



10 Conclusions and outlook 

In the present paper we have reviewed our recent results on the darkon-fluid model which is, 
as outlined in section 9, a first building block for a new general covariant theory describing the 
dark sector of the Universe. The model involves new physics by using nonrelativistic massless 
particles with a nonstandard coupling to the gravitational field. The model, which contains 
no free parameters in its Lagrangian, predicts qualitatively correct values of the late time 
cosmic acceleration as well as the flat behaviour of galactic rotation curves. These successes 
encourage us to continue our work. We have to work out the astrophysical implications of 
the covariant theory and have to compare them with the predictions of the present model and 
with observational data. 

So the main open topics of research are: 

• Derive the EOMs of the covariant theory. 

• Restrict the EOMs to the cosmological regime and compare their predictions with the 
present model. Perform a best fit to the Hubble parameter resp. distance moduli data. 
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• Enlarge the covariant model by adding baryonic matter. 

• Look for the solutions of the spherically symmetric but time dependent EOMs which 
attain for large distances their cosmological limits. 

• Apply the results to 

— modeling of halos, 

— solar system tests. 

• Quantization of the Hamiltonian dynamics for the darkon fluid as a whole or for the 
cosmological dynamics. 

• Answer the question whether the present model is the nonrelativistic limit of the covariant 
theory in a strict sense. 

• Discuss higher order derivative terms within the covariant theory. 

These problems are currently under investigation. We hope to be able to present our results 
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